INVERSE PROBLEMS FOR PERIODIC 
GENERALIZED JACOBI MATRICES 



MAXIM DEREVYAGIN 



Abstract. Some inverse problems for semi-infinite periodic generalized Ja- 
cobi matrices are considered. In particular, a generalization of the Abel cri- 
terion is presented. The approach is based on the fact that the solvability of 
the Pell-Abel equation is equivalent to the existence of a certainly normalized 
J- unitary 2 X 2-matrix polynomial (the monodromy matrix). 



1. Introduction 



In 1826, Abel proved that the square root vi? of the polynomial R of degree 2n 
without multiple zeros can be expanded into a periodic continued fraction, that is 

1 11 11 



(1.1) y/R-U 



Po + j- 

Pi + - 



l p + Pl + ••• + Ps-l + PO 



where U and po, • ■ • , p s -i are polynomials, if and only if the Pell-Abel equation 
(1.2) X 2 - RY 2 = 1 

has the polynomial solutions X and Y (the Abel criterion). Actually, the basic idea 
of Abel was to find out if there exists a polynomial p such that the integral 

VW) 

can be expressed in terms of the elementary functions. It turned out that such 
a polynomial p exists if and only if the Pell-Abel equation (ll.2[) is solvable in 
polynomials. The proof of these results and much more modern information on 
this topic from the algebraic point of view can be found in |15j . 

Another appearance of the Pell- Abel equation is intimately related to extremal 
polynomials. For example, the Chebyshev polynomials Tj and Uj of the first and 
second kind satisfy the following relation 

T 2 (\)-(\ 2 -l)U 2 (\) = l. 

This observation can be generalized to the case of the polynomials that least deviate 
from zero on several intervals [16] . Furthermore, by setting 

n— 1 

R(X) = (X 2 ~l)Y[(X-a J )(X-(3 3 ), 

j=i 
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where — 1 < ai < fi\ < ■ ■ ■ < a n -i < fin-i < 1, we have that the equation fj 1 . 2[) is 
solvable in polynomials if and only if the set 

n-i 

[-1,1] \ |J(«i.A0 

coincides with the spectrum of a bi-infinite periodic classical Jacobi matrix [16] . 
In turn, such a spectrum gives rise to periodic solutions of the Toda lattice (for 
example, see [T71 Chapter 12]). 

On the other hand, we have seen semi-infinite classical Jacobi matrices at the 
beginning although it was implicitly. Indeed, in some special cases the continued 
fraction (II. ip can be a J- fraction generating a classical Jacobi matrix pQ. How- 
ever, J-fractions do not cover all the possible cases of periodic fractions of the 
form Ijl.ip . That is why it is more natural to consider semi- infinite generalized Ja- 
cobi matrices associated with continued fractions of the form (|1.1[) . Such matrices 
were introduced in [8] and the theory was further developed in [7J[9]. In particular, 
direct spectral problems for the periodic generalized Jacobi matrices were consid- 
ered in [7]. The goal of the present paper is to solve some inverse problems for 
semi-infinite generalized Jacobi matrices associated with periodic continued frac- 
tions. It is done in the following way. First, a one-to-one correspondence between 
certainly normalized 2 x 2-matrix polynomials {the monodromy matrices) and the 
periodic generalized Jacobi matrices in question is established. Then we show that 
the solvability of the Pell- Abel equation is equivalent to the existence of a nor- 
malized J-unitary 2 x 2-matrix polynomial, that is the monodromy matrix of the 
underlying periodic generalized Jacobi matrix. The latter statement generalizes the 
Abel criterion. Namely, we give necessary and sufficient conditions for functions of 
the form 

V(X) 

to be the m-functions of semi-infinite periodic generalized Jacobi matrices or, equiv- 
alently, to admit the periodic continued fraction expansions. 

Finally, it should be mentioned that some inverse spectral problems for finite 
generalized Jacobi matrices were studied in [6]. 

2. Preliminaries 

2.1. P-fractions. Let ip be a nonrational function holomorphic at a neighborhood 
of infinity and having the property 

So, ip has the following representation in a neighborhood of infinity 

oo 

(2-1) *W = -£$r. 

where the moments Sj are real. A number rij £ N is called a normal index of the 
sequence s := {sj}°°^ if det(si + fc)"^ = Jj =^ 0. Since tp is not rational, there exists an 
infinite number of normal indices of s (see [TUl Section 16.10.2]). Also, without loss 
of generality we will always assume that the sequence s is normalized, i.e. the first 
nonvanishing moment s ni _i has modulus 1. 
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(2.2) 



As is known [14] , the series (|2.1|) leads to the following infinite continued fraction 

l 



Po(A) - pi(X) - ■■■ - pj(X) 



where Sj = ±1, bj > andpj(X) = X kj +Pk j -\X kj * + ' ' '+Pi A+Po are rea l mon i c 
polynomials of degree fcj (see also [5], [5]). Note, that ; = k + k\ + ■ ■ ■ + kj-\. 

The continued fraction (|2.2p is called a P-fraction. Actually, the P-fraction can 
be considered as a sequence of the linear- fractional transformations [121 Section 5.2] 

Ti(w):= 



p i (A) + £ i 6ja; 
having the following matrix representation 

(2-3) W j (X)=(° eS)' J' GZ +- 

The superposition To o T\ o • • • o Tj of the linear-fractional transformations corre- 
sponds to the product of the matrices W; (A) 

(2.4) W [0J] (A) - K^(A))? fc=1 := W„(A)Wi(A) . . . Wj(X). 

It is well known that the entries of W[oj] can be expressed in terms of denom- 
inators and numerators of the convergents to the P-fraction (|2 . 2[) . To give these 
formulas explicitly define the polynomials Pj +1 (A), Qj + i(X) by the equalities 

The relation W [0j] (A) = Wp^-i] (A)VV,- (A) yields 

« M (A> (J)- W[oj-i](A) ( 4 ) . ( ) . iER 

So, the matrix Wroji (A) has the form 

Further, the equality 

( ~p£$) ) = w ^(xm(X) ( J ) = (A) ( -g } ) , j e N, 

shows that the polynomials Pj(A), Qj(X) are solutions of the difference equation 

Ej-iSjbj-iUj-i —pj{X)u 3 + bju ]+ i — j e N, 
obeying the initial conditions 

Po(A) = l, ft(A) = ^, 

Q (A) = 0, Q 1 (A) = f. 

According to ()2.5|) . the (j + l)-th convergent of the P-fraction (12. 2|) is equal to 
fi := To oTr o • • • o T, (0) = -Q j+1 (A)/P i+1 (A). 
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The relations flU}, ((231), and im ply the following relation [9] 
(2.6) eMQi+i (A)Pf(A) - Q,(A)P J+1 (A)) = 1, j e Z+. 

2.2. Generalized Jacobi matrices. Let p(A) = PnA" + • • • + P1A+P0 be a monic 
scalar real polynomial of degree n, i.e. p n = 1. Let us associate with the polynomial 
p its symmetrizator B p and let the companion matrix C p be given by 

/0 ... -po \ 



/Pi 



\Pn 



"Pi 



\0 1 -Pn-l/ 

As is known, dct(A — C p ) = p(A) and the matrices E p and C p are related by 
(2.7) 



CpEp — EpCp 



Now we are in a position to associate a tridiagonal block matrix with (12.2H . 

Definition 2.1 ([8], |13|). Suppose we are given a P-fraction of the form (|2.2 
Lei 6e rea/ scalar monic polynomials of degree kj 



Pj (\) = \"i + pW)_ A *i- A + • • • + p^A + p^, 



0')' 



and let Sj — ±1, feo > 0, j € Z_| 



TTie tridiagonal block matrix 

Bn 



(2.8) 



where A« 



(2.9) Bj 



Ai Si 
Bi A 2 

Vo ••• -7 

C Pj and fcj+i x fcj matrices Bj and kj x matrices Bj are given by 



'0 



'0 



-B, 



v0 











to// 6e called a generalized Jacobi matrix associated with the P-fraction (|2.2p . 

Remark 2.2. The papers [5], [5], and [13] are only concerned with the case of 
generalized Jacobi matrices which are finite rank perturbations of classical Jacobi 
matrices. 

Now, introducing the following shortened matrices 

(A Q B \ ( A x Bi 



(2.10) H [0Ji 



B Q Ai 



B 



\ 



B\ A 2 



B 



Bj-l A, J \ Bj ; ^ / 

where j G Z + U {oo}, one can obtain the connection between the polynomials Pj, 
Qj and the shortened Jacobi matrices H[o,j], H[i,j]i respectively [8]: 



Pj(X) = (b .../' J -i)" 1 det(A-B [0j _ 1] ), 
?j(A) = £ (b . • . ^-i) -1 det(A - H [hj _ 1} ). 
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Furthermore, the following statement holds true. 

Proposition 2.3 (MM)- Let 3 € N - Then 

i) The polynomials Pj and Pj+\ have no common zeros. 

ii) The polynomials Qj and Qj+\ have no common zeros. 
hi) The polynomials Pj and Qj have no common zeros. 

In what follows we are only interested in periodic generalized Jacobi matrices. 

Definition 2.4. Let s € N. A generalized Jacobi matrix satisfying the properties 

Aj S+ k = A k , Bj S+ k = Bk, £j S+ k = £k, i e k e {0, . . . , s - 1} 

will be called an s-periodic generalized Jacobi matrix. 

Let £? Q x denote the Hilbert space of complex square summable sequences 
(wq, Wi, . . . ) equipped with the usual inner product. Evidently, any s-periodic gen- 
eralized Jacobi matrix generates a bounded linear operator in £? +oo ^ . To see some 
more properties of periodic generalized Jacobi matrices let us define the symmetric 
matrix G by the equality 

(2.11) G = diag(Go,Gi,...) ) Gj = SjE^ 1 , jeZ+. 

Clearly, the operator G defined on £^ Q ^ is bounded and self-adjoint. Moreover, 
G _1 is a bounded linear operator in £? 

Let i0[o,oo) be a space of elements of £? ^ provided with the following indefinite 
inner product 

(2-12) [x,y] = (Gx,y)^ oo) , x,yeif 0iOo) . 

Let us recall H] that a pair (Sj, [•,•]) consisting of a Hilbert space and a 
sesquilinear form [•, •] on fj x is called a space with indefinite inner product. A 
space with indefinite metric (Sj, [•,•]) is called a Krein space if the indefinite scalar 
product [•, •] can be represented as follows 

[x,y] = (Jx,y)sj x,y€fi, 

where the linear operator J satisfies the following conditions 

J= J- 1 = J*. 

The operator J is called the fundamental symmetry. So, one can see that the space 
#[0,00 ) is the Krein space with the fundamental symmetry J = signG (see [1] for 
details). Moreover, the property (|2.7I) implies the following. 

Proposition 2.5 ( 7 , 8 ). The s-periodic generalized Jacobi matrix defines a bounded 
self-adjoint operator H in the Krein space #[o.oo); that is, 

[Hx,y] = [x,Hy] x, y e #[ 0: oo)- 

The main tool for the spectral analysis of periodic generalized Jacobi operators 
is the following matrix 

-£ s _i& s _i(2 s _i(A) -Q S (A) N 
£ s _i& s _ 1 P s _ 1 (A) P S (A) 

The matrix T is called the monodromy matrix. In fact, using the Floquet theory [7] 
we can get the description of the spectrum of the matrix H . Indeed, let wi = wi (A) 



(2.13) T(A):=W [0>s -i ] (A) = f" 
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and W2 — W2(X) be the roots of the characteristic equation det(T(A) — w) = 0. 
Introducing the following notations 

E:={XeC:\w 1 {X)\=\w 2 (X)\}, 

E p := {A e C : P a _i(A) = 0, |6 s _iQ s _ x (A)| > |P S (A)|}. 

we can give the description of spectra of periodic generalized Jacobi operators. 

Theorem 2.6 ( 7 ). The spectrum o~(H) of the s-periodic generalized Jacobi oper- 
ator H has the form 

<t(H) = EUE p , <j p {H) = E p , 
where o~ p (H) denotes the point spectrum of H, i.e. eigenvalues. 

It should be mentioned that for the classical Jacobi matrices this result was 
obtained in 

Remark 2.7 ( 7 ). One can also easily get another description of E 

E = {XeC: (P S (A) - e»_ 1 6»Q a - 1 (A)) e [-2, 2]}. 

2.3. m-functions of periodic generalized Jacobi matrices. Recall that the 
m-functions are one of the central tools in studying linear difference operators. We 
start with the definition of the m-function. 

Definition 2.8. Let H be a bounded generalized Jacobi matrix. The function m 
defined by 

(2.14) m(X) = [(H-X)- l e,e], e = (1, 0, 0, . . . ) T 

is called the m-function (or, the Weyl function) of the operator H . 

Next, by using the Frobenius formula, one can see that the m-function m and 
the m-function 

m-[i,oo) (A) = [(-H[i,oo) - A) _1 e, e] = (G[ 1)00 ) (-%,«,) - A) _1 e, e) 

of -Hhoo) (f° r the definition of ff[i,oo) ancl ^fi,oo) see (|2.10p ) are related by the 
equality [8] 



(2.15) m(A) = - " --, \\\>\\H\\>\\H [hoo) \ 



-£q 

p (X) +£o6o TO [i,oo)(A)' 

The standard technique of dealing with periodic continued fractions enables us 
to calculate the m-function explicitly. 



Proposition 2.9. Let H be an s-periodic generalized Jacobi matrix. Then its m- 

function has the following form 

(2-16) 

... -(P S (A) + e._i6._iQ._ 1 (A)) + y/(PsW - e s _i& s _ig s _i(A)) 2 - 4 

m(X) = — , 

V j 2e s _ 1 & s _ 1 P s _ 1 (A) 

where the cut is E and the branch is chosen in such a way that m(A) — > as X — > oo. 
Proof. Due to (|2.15p and the s-periodicity of H, we have that 

7 , 2 



m(A) = — 



e £o£i^o £ s -i£ s - 2 b 2 s 



p {X) - pi(X) - ■■■ - Ps-i(X) + e s -i^_ 1 m(A) 



INVERSE PROBLEMS 



7 



According to (|2.4|) and (|2.13|) the latter relation takes the form 

, 9iri m _ £ s -ib s -iQa_i(A)m(A) + Q S (A) 

1 ' [) ~ e s ^ 1 P s . 1 (X)m(X) + P s (X)- 

Next, (|2.17[) can be rewritten as follows 

(2.18) £ s _ 1 6 s _ 1 P s _ 1 (A)m 2 (A) + (P. (A) + e,_i& a _iQ._i(A))m(A) + Q S (X) = 0. 

Now, (I2.16[) follows from (|2.18p and (|2.6p . The choice of the cut and the branch is 
due to Definition 12.81 and Theorem 12.61 □ 



Theorem 2.10 (cf. [7 ). The m-function of the periodic generalized Jacobi matrix 
H determines H uniquely. 

Proof. In fact, consecutive applications of the relation (|2.15[) lead to the continued 
fraction (|2.2I) . which determines H uniquely. □ 

3. Inverse problems 

3.1. The monodromy matrices. Note that knowing the monodromy matrix T 
gives the complete information about the spectrum of the corresponding s-periodic 
generalized Jacobi matrix H. Moreover, it is clear that one can reconstruct the 
s-periodic generalized Jacobi matrix H by its monodromy matrix T = W[o,«— i] 
since by the construction we have that 

Q S (X) e e eibl £ s _ 2 £ s _i& 2 _ 2 



P S (X) p (X) — pi(X) — ■■■ - p s -i(X) 

and 6 s _i can be extracted from T directly. Now, one of the most natural questions 
is the following: which matrices can be the monodromy matrices? To answer this 
question let us note first that the monodromy matrix 



rm _ ftn(X) t 12 (X)\ 



satisfies the following properties 

(Tl) t±i, tn, tn, ti2 are polynomials with real coefficients such that 

<n(A)< 22 (A) - t 12 (A)i 2 i(A) = 1, Ael; 

(T2) degt i2 < degt 22 , degi 2i < degt 22 ; 

(T3) the absolute values of the leading coefficients of £21 and £ 22 are equal and 
the leading coefficient of i 22 is positive. 

These properties give us a tip to introduce the following definition. 

Definition 3.1. The 2 x 2-matrix polynomial T is called admissible if it satisfies 
the conditions (Tl)-(TS). 

Remark 3.2. Consider the following 2 x 2-matrix 

It is said that 2x2 matrix polynomial is J-unitary on the real line R if the following 
equality holds true 

T(X)JT*(X) = J, XeR. 
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Clearly, every admissible matrix polynomial is J-unitary. Moreover, every J- 
unitary 2 x 2-matrix polynomial on the real line can be normalized by multipli- 
cations from left and right by constant J-unitary matrices to an admissible matrix 
(for more details see [9j Section 3.7]). 

The answer to the above-stated question is given by the following statement. 

Theorem 3.3 (cf. [9). Every admissible matrix polynomial T is the monodromy 
matrix of a unique periodic generalized Jacobi matrix. 

Proof. Actually, this statement can be proved by re-examining [9J Section 3.7]. 
First, note that the admissible matrix T generates the shortened generalized Jacobi 
matrix H[ jS _i]. Indeed, let us consider the P-fraction expansion 

ti2(A) _ £p £q£i6q e s -2e s -ib 2 s _ 2 
t 22 (A) ~ po(A) - Pi (A) - ••• - p s _x(A) ' 
which, according to the definition, generates the shortened generalized Jacobi ma- 
trix -ffro.s-i]- So, one can recover the sequences Po, Pi, . . . , P._i, 6 s -iP s and Qq, 
Qx, ■ • • , Qs-i, 6 s -iQ s such that tnjt-n = —Q s /P s . Choosing 6 s _i > to have the 
equalities t±2 — ~Qs and t22 = Ps we conclude that we know the matrix 

-e s _i&._iQ,_i (A) -Q S (\T 
e._i6._iP._i(A) P.(A) 

and the number 6._i. Also, the periodicity implies that £o = s s and, so, one can 
determine = £._i£.6._i. Thus, the admissible matrix T gives rise to the 

following s-periodic generalized Jacobi matrix 

/ff [0lS -i] B.-i 0\ 

B s ~l ff[o,._l] Ps-l 

B s -i P[o, s -i] 

V o 

Finally, it follows from (T1)-(T3), Proposition 12.31 and (|2.6p that the monodromy 
matrix Wro,._i] of H coincides with T (for more details see the proof of [9[ Theo- 
rem 3.23]).' □ 

Remark 3.4. Due to Remark 13.21 and formula (|2.4[) . Theorem [33] can be applied for 
getting factorization results for J-unitary matrix polynomials [2] , [9] . 

3.2. The Abel criterion. Formula (12.16[) leads to the conclusion that the m- 
function of the periodic generalized Jacobi matrix has the following algebraic form 



W[o,.- 



-i](A)=(" 



H 



(3-D ,(A, . ^ m Um , 

where R, U, and V are polynomials such that 

(3.2) deg R = 2n, deg U = n, deg V < deg U. 

Besides, the m- function has the property that ip(\) — > as A — > oo. 

Clearly, the function if can be also represented by means of some other choice 
of polynomials R, U, and V . So, without loss of generality we always assume that 
the polynomials R, U, and V are chosen in such a way that the polynomial V 
has the minimal degree between all such polynomials in the representation (|3.1[) of 
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the given function (p. Now the natural question is to ask under what conditions 
on R, U, and V a function of the form (|3.1|) can be the m-function of a periodic 
generalized Jacobi matrix? One can find the answer in the following statement. 

Theorem 3.5. Let R, U , and V be polynomials subject to (|3.2p . Let us also suppose 
that one of the branches of the function ip — (y/R — U)/V possesses the property 
p(X) — > as A — > oo. Then there is a cut such that the corresponding branch of ip 
is the m-function of a periodic generalized Jacobi matrix if and only if there exist 
real polynomials X , Y , Z satisfying the following relations 

(3.3) X 2 - RY 2 = 1, 

(3.4) {U 2 - R)Y = VZ. 

Proof. The "if" part. Let <p = (y/R — U)/V be the m-function of a periodic gener- 
alized Jacobi matrix. Then ip satisfies the following equation 

(3.5) e._i6._iP._i(A)^ 2 (A) + (P.(A) + e,_i&,_iQ a _i(A)MA) + Q S (A) = 0. 
Substituting tp = (y/R - U)/V to ([33) one gets 

, s VR(-2e s -ibs-iPs-iU + P s V + e s -ib s - 1 Q s - 1 v) + 

(3.6) V J 

+e s _ 1 6 s _ 1 P s _i(i? 2 + U 2 ) - (P s + e,_ift,_iQ 5 _i)l/V + Q S V 2 = 0. 

Observe that the equality ay/~R -1-/3 = 0, where a and j3 are some polynomials, 
implies a = and /3 = 0. So, it follows from (|3.6p that 

(3.7) P S V = 2e._i6._iP._il7 - e a _ 1 6._ 1 Q a _ 1 V ) 

(3.8) Q S V 2 = {Ps+es-ibs^Qs^UV -e s -ibs-iPs-i{R 2 + U 2 ). 
Combining (|3.7I) and (|3.8I) gives 

(3.9) Q.y 2 = e._i6._iP._i(C/ 2 -P). 

It is easy to see from (I2.16|) as well as from (I2.18[) that P s _i = P._iV . Taking this 
observation into account one can rewrite (I3.7[) and (I3.9[) as follows 

(3.10) P s = 2e._i6._iP._iC7-e._i6._iQ._x, 

(3.11) Q S V = e._i6._iP._i(C/ 2 -P). 

Now, by setting Z = Q s and Y = e._i6._iPs_i the relation (|3. 1 1|) becomes (|3.4p . 
In order to get (|3.3[) . let us multiply formula (|2.6[) by y 

e._i6._i (FQ s P s -i - FQ s -iP s ) = 1. 

Substituting (|3.10[) and p. lip to the latter equality we arrive at 

e._i6._i (es-ibs-iP^ViU 2 - R) - e._i6._i7(2t/P._i - Q._i)Q._i) = V. 
Further calculations lead to the relation 

bl-i (Ps-i(U 2 -R)- 2UP s -iQs-i + Ql-x) = I, 
which can be rewritten as follows 

(3.12) 6 2 _ 1 (P s _iC/ - Q s ^) 2 - 6 2 _ 1 PP. 2 _ 1 = 1. 

Setting X = e s -ib s -i(P s -iU — Q.-i) and recalling that Y — £ s _ife s _iP s _i one 
sees that p.!2p is exactly p.3p . 
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The "only if" part. Suppose that the system ()3.3|) . (|3.4j) is satisfied for some 
polynomials X, Y, and Z. Then the "if" part gives us a hint to introduce the 
following polynomials 

tn = X - YU, iia = YV, iai = -Z, t 22 = A + Ff/, 

where the signs of the leading coefficients of X, Y, and Z are taken to have the 
leading coefficient i 22 positive and degin < n. Next it is easy to check 

det T = iutaa - iiaiai = X 2 - Y 2 U 2 + ZYV = 

= X 2 - (RY 2 - ZYV) + ZYV = X 2 - RY 2 = 1. 

It remains to note that 

degt 2 i = degZ = deg{U 2 - R)Y/V < dcgX = degi 22 

in order to see that the matrix T = (ty)? - =1 is admissible. The rest follows from 
Theorem ELl □ 

Remark 3.6. If V = 1 then Theorem 13.51 reduces to the classical Abel criterion of 
the periodicity of continued fractions representing y/~R. 

Remark 3.7. Every solvable system (|3.3[> . (|3.4|l leads to an admissible matrix poly- 
nomial. The converse is also true. In turn, every admissible matrix polynomial can 
be constructed by $273]), (|2^4| . 
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